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I Abstract 

The currency carry trade is the investment strategy that involves selling low interest rate currencies 
^yy I in order to purchase higher interest rate currencies, thus profiting from the interest rate differentials. 

(Ih . This is a well known financial puzzle to explain, since assuming foreign exchange risk is uninhibited 



and the markets have rational risk-neutral investors, then one would not expect profits from such 
strategies. That is uncovered interest rate parity (UIP), the parity condition in which exposure to 



' foreign exchange risk, with unanticipated changes in exchange rates, should result in an outcome that 

>'. 

, changes in the exchange rate should offset the potential to profit from such interest rate differentials. 

. The two primary assumptions required for interest rate parity are related to capital mobility and 

' perfect substitutability of domestic and foreign assets. Given foreign exchange market equilibrium, 

I the interest rate parity condition implies that the expected return on domestic assets will equal the 

exchange rate-adjusted expected return on foreign currency assets. 

However, it has been shown empirically, that investors can actually earn arbitrage profits by bor- 
rowing in a country with a lower interest rate, exchanging for foreign currency, and investing in a 
foreign country with a higher interest rate, whilst allowing for any losses (or gains) from exchanging 
back to their domestic currency at maturity. Therefore trading strategies that aim to exploit the inter- 
est rate differentials can be profitable on average. The intention of this paper is therefore to reinterpret 
the currency carry trade puzzle in light of heavy tailed marginal models coupled with multivariate tail 
dependence features in the analysis of the risk-reward for the currency portfolios with high interest 
rate differentials and low interest rate differentials. To achieve this analysis of the multivariate extreme 
tail dependence we develop several parameteric models and perform detailed model comparison. 

Keywords: Currency carry trade. Multivariate tail dependence, Forward premium puzzle. Mixture 
models, Generalized Archimedean copula and Extreme value copula. 
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1 Introduction 



In the finance literature there has been a growing tendency to study relationships between different 
financial assets. For example, futures and spot, forward and spot relationships, relationships between 
interest rate products of different maturity, or relationships between individual stocks and their associated 
options etc. Primarily the emphasis of such studies has focused on the role of extreme events in the 
valuation and understanding of such relationships. Until recently the majority of such financial studies 
have undertaken such analysis from the perspective of modelling marginal behaviours through the analysis 



of extreme value model 
length review in 



ower law relationships and heavy tailed models, see for example a book 



RachevI ()2003l ). In this paper we propose to utilize the ideas from such studies for the 
marginal model analysis whilst also developing multivariate joint tail dependence modelling components. 
The aim will be to utilise such analysis to revisit and reinvestigate, in light of such multivariate tail 
dependence models, a well known puzzle in finance which is known as the currency carry trade or 
violation of the uncovered interest rate parity condition (UIP). 

The UIP parity condition, discussed in detail below, provides a unique relationship that has implica- 
tions for the forecasting of future exchange rates: the forward exchange rate and the future spot exchange 
rate. In this paper we aim to explore aspects of trading strategies that have been shown to violate the 
assumptions of UIP and are able to obtain profits from portfolios constructed around interest rate dif- 
ferentials. In particular we aim to explore to what extent one can attribute these excess average returns 
with regard to compensation for exposure to tail risk, i.e. dramatic depreciations in the value of the high 
interest rate currency in times of high FX volatility. 

We postulate that such analyses should also benefit from consideration not only of the marginal be- 
haviours of the processes under study, in this case the exchange rates for currencies in a portfolio, but 
also a rigorous analysis of the joint dependence features of such relationships. Therefore we investigate 
such joint relationships in light of the UIP condition. We argue that the analysis of the carry trade is 
better informed when the model analysis is performed jointly by modelling the multivariate behaviour of 
the marginal processes accounting for potential multivariate extremes, and multivariate tail dependence 
features, whilst still incorporating heavy-tailed relationships studied in marginal processes. To achieve 
the accurate modelling of the potentially complex multivariate dependence features the statistical models 
considered should be sufficiently flexible to accommodate a variety of potential tail dependence relation- 
ships. In this regard we consider two flexible families of mixture copula models comprised of members of 
the Archimedean copula family and their outer-power transform extensions. 

We propose to utilise such flexible multivariate statistical models to reinvestigate from a new perspec- 
tive involving copula dependence features the currency carry trade portfolio behaviour. Such trading 
strategies are a popular approach which involves constructing highly leveraged portfolios which involve 
selling low interest rate currencies in order to buy higher interest rate currencies, thus profiting from the 
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interest rate differentials. When such portfolios are highly 



leveraged th is can result in sizable profits . 



Famal (jlDSJ) or iHansen and Hodrickl (jl98d ) 



The presence of such arbitrage opportunities pointed out by 
violates the fundamental relationship of uncovered interest rate parity (UIP). The UIP refers to the parity 
condition in which exposure to foreign exchange risk, with unanticipated changes in exchange rates, is 
uninhibited and therefore if one assumes rational risk-neutral investors, then changes in the exchange 
rate should offset the potential to profit from such interest rate differentials. The existence of a working 
UIP relationship is related to two primary assumptions which are the capital mobility and the perfect 
substitutability of domestic and foreign assets. When UIP holds, then given foreign exchange market 
equilibrium, the interest rate parity condition implies that the expected return on domestic assets will 
equal the exchange rate-adjusted expected return on foreign currency assets. Therefore, no arbitrage 
opportunities should arise in practice; however such opportunities are routinely observed and exploited 
by large volume trading strategies. 



In this paper we aim to explore aspects of simple trading strategies that have been shown to exploit 
this violation of the UIP assumptions and are therefore able to obtain profits from portfolios constructed 
around interest rate differentials. In particular we investigate whether such strategies, though seemingly 
profitable are actually exposed to significant risk that can only be observed in the presence of suitable 
models for the joint extreme tail dependence in such portfolios. Therefore, we model the extent of these 
excess average returns and discuss the potential risk-reward relationships when one adequately accounts 
for the modelling of the extreme joint tail-dependence relationships between such currency exchange rates 
used to construct the portfolios. We demonstrate the sustained presence of multivariate extreme tail de- 
pendence features between the currencies typically used to construct the high-interest rate portfolios and 
the low interest rate portfolios in practice. In addition we observe evidence for different structural forms 
of the observed tail dependence in periods of high market volatility compared to periods of low volatility. 
This is observed when comparing multivariate tail dependence features of a portfolio with high interest 
rate differentials against the tail dependence features of a portfolio with low interest rate differentials. 



This phenomenon will have implications on the assessment of the profitability of such portfolios when 
one takes into account suitable consideration of the compensation required for exposure to tail risk, i.e. 
joint dramatic depreciations in the value of the high interest rate currencies in times of high FX volatility. 
Therefore, our analysis should provide interesting insight into the seemingly complex understanding of 
how on average arbitrage profits are regularly made on highly leveraged portfolios, however the presence 
of multivariate tail dependence features in the assets of the carry portfolio will change the manner in 
which one should assess the risk associated with such highly leveraged profits. This may help to ex- 
plain the violation of the UIP condition when such risk analysis is appropriately considered for the joint 
behaviour. It is our conjecture that such profits and the risk metrics considered when assessing the risk- 
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return relationships of such trading strategies should be therefore altered to account for tail dependence 
features in the portfolios constructed from interest rate differentials. We would therefore argue that 
failure to take into consideration the potential for significant losses associated with a non-trivial joint 
probability of large adverse movements in exchange rates could wipe out all potential average profits 
from a single event should it occur. 

To summarise, in this paper we aim to study the multivariate relationships with respect to upper and 
lower joint tail-dependence within currency carry trade portfolios which are constructed from portfolios 
based on the actual interest rate differentials between currencies. The paper is structured as follows. In 
the second section, we review the concept of UIP and the related literature. The third section presents 
details of the portfolio construction as well as a detailed description of the data utilised in the analysis. 
The fourth section is devoted to the description of the copula models and their applications in finance. 
We also describe the model we propose to measure the tail dependence within the high interest rate and 
low interest rates portfolios. Section five presents estimation of the multivariate distributions for the 
mixture copula and marginal heavy tailed models via the Inference on the Margins approach. Finally 
section six and seven provide results and comments about the relationships between tail dependence and 
the average performance of the different currency portfolios and conclusions. 



2 Currency Carry Trade 



Numerous empirical studies (jFama 



19841 : 



Hansen and Hodrick 



1980 



Engel, 



1993; 



Lustig and Verdelhan 



20071 ) have previously demonstrated, that investors can actually earn arbitrage profits by borrowing in a 
country with a lower interest rate, exchanging for foreign currency, and investing in a foreign country with 
a higher interest rate, whilst allowing for any losses (or gains) from exchanging back to their domestic 
currency at maturity. Therefore, trading strategies that aim to exploit the interest rate differentials can 
be profitable on average. This is notably the case for the currency carry trade which is thus the simple 
investment strategy of selling a low interest rate currency forward and then buying a high interest rate 
currency forward. The idea is that the interest rate returns will outweigh any potential adverse moves 
in the exchange rate. Historically the Japanese Yen and Swiss Franc have been used as "funding curren- 
cies" , since they have maintained very low interest rates for a long period. The currencies of developing 
nations, such as the South African Rand and Brazilian Real have been typically used as "investment 
currencies" . Whilst this sounds like an easy money making strategy there is of course a downside risk. 
This risk comes in the form of currency crashes in periods of global FX volatility and liquidity shortages. 
A prime example of this is the sharp Yen carry trade reversal in 2007. 
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In this section we first present more formally the forward premium puzzle and the aforementioned 
currency carry trade strategy that can be implemented to benefit from this violation of the UIP relation. 
Then, we detail some of the justifications that have been recurrently proposed in the literature before 
proposing another explanation based on the analysis of the tail-dependence between the high interest 
rates and the low interest rates currencies. 



2.1 The Forward Premium Puzzle 



'his 



pheno mena introdu ced initially by 



Hansen and Hodrickl (|l98Cll ) , [Hansen and Hodrickl (jl983l ) , 



Fama 



(| 19841 ) and lEngell (| 19841 ) is directly linked to the arbitrage relation existing between the spot and the 
forward prices of a given currency, namely the Covered Interest Parity. This relation states that the price 
of a forward rate can be expressed according to the relationship: 

irt-r{){T-t) g 



FT = 



(2.1) 



where Ff and St denote respectively the forward and the spot prices at time t. While rt and r/ represent 
the local risk free ratqj and the foreign risk free rate. We denote by T the maturity of the forward 
contract considered. It is worth emphasizing that under the absence of an arbitrage hypothesis, this 
relation is directly resulting from the replication of the forward contract payoff using a self financed 
strategy. Moreov er, it has been demonstrated ernpiric ally the validity of this arbitrage relation in the 



currency market (jJuhl et al 



20061 : 



Akram et al. 



20081 ) when we consider daily data. Then, if we take 



the logarithm of expression (|2.1|) we thus obtain the following relation: 

fl - ST = {n - r{)T - {sT - St) 
= in - r{)T - AsT 



(2.2) 



where r is the time to maturity from calendar day t given by (T — t) while and st denote respectively 
the log-values of the forward and the spot prices. 

This expression notably represents the excess return received by an investor who has invested one 
unit of local currenc}0 in a forward contract of maturity T at time t and has held this position until 
expiration to convert it back to his reference currency with an exchange rate equal to St- We can notice 
that at the trade settlement the profit or the loss in local currency for this investor equals the differential 
of interest rates on a prorata temporis basis, plus the differential of exchange rate betwe en t and T. 



Finally, if we assume the forward price is a martingale under the risk neutral probability (jMusiela and Rutkowsk: 



2011 



) then its value equals: 



EQlSrlJ't] = F^ 



(2.3) 



^We mean by local risk free rate the interest rate prevailing in the reference country which would be for instance the 

dollar for an American investor. 

^Which corresponds in our case to the US dollar since we consider the position of an American investor. 
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where J-t is the filtration associated to the stochastic process St- Replacing the expression (|2.3p in the 
relation ()2.ip leads to the formula: 

= EL = e(n-rf)(T-t) ^2.4) 
St 

Which leads then to the UIP hypothesis since according to this expression, under the risk neutral 
probability the expected variation of the exchange rate St should equal the differential of interest rate 
between the two countries. Thus if an investor buys a forward contract the profit or the loss resulting at 
the contract maturity should accordingly be equal to zero as the exchange rate at maturity St should be 
equal to Ste^^^~^'t)(.'^-'^) which was the price paid initially for the forward contract under the hypothesis 
of absence of arbitrage opportunities. 



St 



St 



2.2 A Review of the Literature 



As discussed above, if the UIP relationship held then there should indeed not be any yield difference 
between a risk-free investment in a reference currency and a risk-free investment in another currency 
after converting it back to the reference currency. Accordingly, the depreciation of a currency rel- 
ative to another should be equal to the risk fre e inter est r ates di f ferent ial between them. However, 



Hansen and Hodrick (198C). 



Lustig et al. 



2011 



Hansen and 



Lustig and Verdelhan 



2007 



lodrick mm and 



Menkhoff et al. 



Fanial (119841 ) among other recent articles 



20121 ). demonstrate that this relation 



is not observed empirically in markets data and that a strategy named "currency carry trade" can even 
benefit from this fiaw. 

Among the justifications of this phenomenon. iFamal (jl984l ) i nitially proposed a time varying risk premia 



WeitzmanI (|2007l ) demonstrates through a 



within the forward rate relative to the associated spot rate. 
Bayesian approach that the uncertainty about the variance of the future growth rates combined with 
a thin-tailed prior distribution would generate the fat-tailed distribution r equired to solve t he fo rward 



Menkhoff et al. 



(120JJ) who 



premium puzzle. This could be compared to the argument retained by 
demonstrate that high interest rate currencies tend to be negatively related to the innovations in global 
FX volatility which is considered as a proxy for unexpected c hanges in the FX market volatility. 



Another hypothesis, proposed by 



Farhi and Gabaixl (j2008l ). consists in justifying this puzzle through 



the inclusion of a mean reverting risk premium. According to their model a risky country, which is more 
sensitive to economic extreme events, represents a high risk of currency depreciation and has thus to 
propose, in order to compensate this risk, a higher interest rate. Then, when the risk premium reverts to 
the mean, their exchange rate appreciates while they still have a high interest rate which thus replicates 
the forward rate premium puzzle. 

The causality relation between the interest rate differential and th e currency shocks can be presented 
the other way around as detailed in lBrunnermeier and PedersenI ((20091). In this article, the authors indeed 
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assume that the currency carry trade mechanically attracts investors and more specifically speculators 
who accordingly increase the probability of a market crash. Tail events among currencies would thus be 
caused by speculators' need to unwind their positions when they get closer to fundin g constraints. 



This recurrent statement of a re lation between tail events and forward rate premium (|Farhi and Gabaix 



20081 : 



Brunnermeier et al 



2008 ) has led us to propose in this article a rigorous measure and estimation 
of the tail thickness at the level of the marginal distribution associated to each exchange rate. Moreover, 
the question of the link between the currencies marginal distributions and the associated interest rates 
differential leads us to consider more globally the joint dependence structures between the individual 
marginal cdf tails with respect to their respective interest rate differential. 



2.3 Tail Dependences and Forward Premium Puzzle 



The approach we adopt in this paper is a statistical framework with a high degree of sophistication, 



ogous in nature to t 



re ideas considered 



Mehra and Prescott 



(jl985l ) in the late 



however its fundamental reasoning and justification is indeed anaL 
when investigating the "equity risk premium puzzle" coined by 
80's. The equity risk premium puzzle effectively refers to the fact that demand for government bonds 
which have lower returns than stocks still exists and generally remains high. This poses a puzzle for 
economists to explain why the magnitude of the disparity between the returns on each of these asset 
classes, stocks versus bonds, known as the equity risk premium, is so great and therefor e implies an 



implausibly high level of investor risk aversion. In the seminal paper writte n 



proposes to explain the "equity risk premium puzzle" (jMehra and Prescott 



Rietd (j 19881 ). the author 



19851 ) by taking into consid- 



eration the low but still significant probability of a joint catastrophic event. 



Analogously in this manuscript, we are proposing to explore the highly leveraged arbitrage oppor- 
tunities in currency carry trades that arise due to violation of the UIP. However, we conjecture that if 
the assessment of the risk associated with such trading strategies was modified to adequately take into 
account the potential for joint catastrophic risk events accounting for the non-trivial probabilities of joint 
adverse movements in currency exchange rates, then such strategies may not seem so profitable relative 
to the risk born by the investor. We propose a rigorous probabilistic model in order to quantify this phe- 
nomenon and potentially detect when liquidity in FX markets may dry up. This probabilistic measure 
of dependence can then be very useful for risk management of such portfolios but also for making more 
tractable the valuation of structured products or other derivatives indexed on this specific strategy. To 
be more specific, the principal contribution of our article is indeed to model the dependences between 
exchange rates using a flexible family of mixture copulas comprised of Archimedean members and their 
outer-power transforms. This probabilistic approach allows us to express joint distribution of the vectors 
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of random variables, in our case vectors of exchange rates log-returns in each basket of currencies, as 
functions of each marginal distribution and the copula function itself. 



Whereas in the literature we mentioned earlier, the tail thickness resulting from the carry trade 
has been either treated individually for each exchange rate or through the measurement of distribution 
moments that may not be adapted to a proper estimation of the tail dependencies. In our paper, we 
propose instead to build on a daily basis a set of portfolios of currencies with regards to the interest rate 
differentials of each currency with the US dollar. Using a mixture of copula functions, we then extract 
a measure of the tail dependencies within each portfolio and finally interpret the results. Among the 
outcomes of our study, we demonstrate that during the crisis periods, the high interest rate currencies 
tend to display very significant upper tail dependence. Accordingly, we thus conclude that the appealing 
high return profile of a carry portfolio is not only compensating the tail thickness of each individual 
component probability distribution but also the fact that they tend to occur simultaneously and lead 
to a portfolio particularly sensitive to the risk of drawdown. Furthermore, we also demonstrate that 
high interest rate currency portfolios can display periods during which the tail dependence gets inverted 
demonstrating when periods of construction of the aforementioned carry positions are being undertaken 
by investors. 



3 Data Description and Portfolios Construction 

In this section we describe the set of data we use for our empirical study and describe the macro-economic 
specificities associated to some of the currencies we considered. Furthermore, we present the method we 
retained in this article to build the portfolios that we combine later on to build a carry trade position. 

3.1 Data Description 

We consider for our empirical analysis a set of 20 currency exchange rates relative to the USD. We indeed 
considered the point of view o f an American investor as th is is generally the hypothesis retained in the 



literature ( Brunnermeier et al 



2003; 



Menkhoff et al 



201^ ). However the same analysis could be carried 



out from any other investor standpoint as the phenomena we will describe does not only depend on a 
specific currency but more on two sets of currencies. These sets of currencies correspond to the high 
interest rate currencies which are used to obtain the highest return (named the "investment currencies" ) 
and the low interest rate currencies which allows for borrowing at a low cost the amount of money 
necessary for this investment (named the "financing currencies"). 

The time series analysed range from 04/01/2000 to 02/01/2013 and comprise the following currencies: 
Euro (EUR), Turkish Lira (TRY), Japanese Yen (JPY), British Pound Sterling (GBP), Australian Dollar 
(AUD), Canadian Dollar (CAD), Norwegian Krone (NOK), Swiss Franc (CHF), Swedish Krona (SEK), 
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Mexican Peso (MXN), Polish Zloty (PLN), Malaysian Ringgit (MYR), Singaporean Dollar (SGD), In- 
dian Rupee (INR), South African Rand (ZAR), New Zealand Dollar (NZD), Thai Baht (THB), South 
Korean Won (KRW), Taiwanese Dollar (TWD), Brazilian Real (BRL). We have been provided, on a 
daily basis, with the settlement prices for each currency exchange rate as well as the simultaneous price 
for the associated 1 month forward contract. Due to differing market closing days, e.g. national holidays, 
there was missing data for a couple of currencies and for a small number of days. For missing prices, the 
previous day's closing prices were retained. 

The reason why we based our analysis upon a constant maturity 1 month forward is twofold. Firstly, 
we do not try in this paper to replicate as realistically as possible a currency carry trade portfolio to see if 
there is a recurrently high average return. The main inconvenience of such analysis comes from the loss of 
data points. As a matter of fact, to build a carry portfolio, the position has to be held until the maturity 
of the forward contract which leads in our case to retain only one point for each month. However, in our 
case we have at our disposal one point per day which makes our analysis of individual tails and their 
interdependences more robust. Secondly, tail behaviour of monthly data is naturally different from the 
tail behaviour of daily data, one reason for this difference is that individual currencies can display a mean 
reversion in the mid-term and thus reduce the amplitude of the movement. 

Among the currencies under scrutiny, some of them have displayed very large variations in the last 
decade mainly for macro-economic reasons. Therefore, we considered it insightful to mention some of 
the most meaningful. The Brazilian Real displays in its time series two important periods of shocks, 
the first in 2001 and the second in 2002. Naturally the first of them was due to the terrorist attacks 
against the world trade center in September. However the Brazilian real has been also impacted by the 
market's concerns of a contagion after the rumours of default of the Argentinian government. The second 
shock on the Brazilian real in 2002 was related to the potential election of the Workers' Party leader Luiz 
Inacio Lula da Silva which prompted concern he might spark a default by overspending to meet promises 
of spurring growth and employment. In 2001, the South African rand slumped 29% after the events of 
September 11 and the market's concern of a global recession and a slump in commodity prices to which 
the South African economy is particularly exposed to. As a third example of a shock in an instrumental 
currency in a carry trade strategy we note the 30% daily loss of the Turkish Lira on the 22nd of February 
2001. This was due to Turkey's decision to abandon the defense of their currency in order to reduce 
the cost of financing lira-denominated debt. It is worth mentioning that we did not remove these data 
points from our time series given that different events may have impacted the other exchange rates at a 
different time but our analysis does not focus only on the tail events associated to a particular currency 
but more on the events impacting simultaneously a set of currencies. 
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3.2 Currency Portfolios Formation 

As described earlier, the currency carry trade results from the differential of interest rates prevailing in 
different countries. By borrowing a certain amount of money in low interest rate countries and investing 
it in high interest rate countries, a recurrent profit can be generated given that the UIP condition is on 
average not satisfied. In order to differentiate the "financing currencies" from the "investment currencies" , 
we start by classifying each currency relative to its differential of risk free rate with the US dollar. We note 
the following basic explanation of the high rates and low rates. In general countries that are considered 
'safe' can borrow at a lower i nterest rate, which may ex plain why historically the US dollar or Swiss 



Franc interest rates were low (jGourinchas and Rev . 



20071 ) while the Turkish Lira rates were historically 



high as this country is not considered as financially secure. 

Moreover we demonstrated in expression (j2.4p that the differential of interest rates between two 
countries can be estimated through the ratio of the forward contract price and the spot price. Accordingly, 
instead of considering the differential of risk free rates between the reference and the foreign countries, we 
build our respective baskets of currencies with respect to the ratio of the forward and the spot prices for 
each currency. On a daily basis we compute this ratio for each currency and then build five portfolios of 
four currencies each. The first portfolio gathers the four currencies with the highest positive differential 
of interest rate with the US dollar. The selected currencies over the period 04/01/2000 to 02/01/2013 
for the high interest rate basket are displayed in Figure [TJ These currencies are thus representing the 
"investment" currencies, through which we invest the money to benefit from the currency carry trade. 
The last portfolio will gather the four currencies with the highest negative differential (or at least the 
lowest differential) of interest rate. As with the high interest rate basket, we also display the low interest 
interest rate currency selections in Figure [2j These currencies are thus representing the "financing" 
currencies, through which we borrow the money to build the currency carry trade. 

Conditionally to this classification we investigate then the joint distribution of each group of currencies 
to understand the impact of the currency carry trade, embodied by the differential of interest rates, on 
currencies returns. In our analysis we concentrate on the high interest rate basket (investment currencies) 
and the low interest rate basket (funding currencies), since typically when implementing a carry trade 
strategy one would go short the low basket and go long the high basket. 



4 Copula Modelling Context 

In this section we provide a brief survey of some relevant key results that are known and typically utilised 
to motivate the development of multivariate models in a constructive manner as would be carried out 
with a copula framework. What is meant by this is that we aim to consider two stages, first the estima- 
tion of suitable heavy tailed marginal models, followed by the estimation of the dependence component 
of the multivariate model. That is we combine model knowledge of suitably fitting marginal models 
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Figure 1: Basket 5 (highest IR) composition. 
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Figure 2: Basket 1 (lowest IR) composition. 
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with hypotheses of potential features of the tail dependence via a copula to obtain a joint multivariate 
distribution. Such an approach has been motivated by the fundamental results developed in lSklad (|l959l ) 
for copula decompositions of multivariate distributions. 

A copula is specified according to the following definition. 
Definition: Copula 

A d- dimensional copula is a multivariate cumulative distribution function C with uniform [0, 1] margins, 
i.e. C : [0,1]"^ ^ [0,1]. 

Now we demonstrate one of the main attractions for practitioners for the use of copula models, via 
an understanding of separation of a multivariate distribution into its marginal distributions and the 
depedence structure between the margins. Sklar's theorem (j4.ip provides the foundation to the study 
of copulae by proving that any multivariate distribution with continuous margins has a unique copula 
representation. 
Sklar's Theorem (1959) 

Consider a d-dimensional cdf F with marginals Fi, . . . , F^. There exists a copula C, s.t. 

F{xi, ...,xd) = C(Fi(xi), . . . , Fd{xd)) (4.1) 

for all Xi € (— oo, oo),i € 1, . . . ,d. Furthermore, if Fi is continuous for all i = 1, . . . ,d then C is unique; 
otherwise C is uniquely determined only on RanFi x • • • x RanF^, where RanFi denotes the range of the 
cdfFi. 

Copula models therefore provide a mechanism to model the marginal behaviour of each currency and 
then separately to focus on developing hypotheses regarding the possible dependence structures between 
the log returns of the forward exchange rates of the currencies in the portfolios which can be tested 
through parameterization of a model via a copula and then a process of model selection. 

4.1 Archimedean Copula Models 

There are many possible copula models that could be considered in the modelling of the multivariate 
dependence features of the currency portfolios. The intention of this analysis was to work with well 
known models which have well understood tail dependence features and are relatively parsimonious with 
regard to the number of parameters specifying the copula. We obtain flexible dependence relationships 
by combining such components into mixture models that allow for a range of flexible tail dependence 
relationships to be studied. 

In particular, we will focus on the well-known class of Archimedean copulae, as defined in ()4.2p . 
since they provide a parsimonious approach that allows for the modelling of various tail dependence 
characteristics. A function -0 is said to generate an Archimedean copula if it satisfies the properties 
below. 
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Definition: Archimedean Generator 

An Archimedean generator is a continuous, decreasing function ip : [0, oo] — t- [0, 1] which satisfies the 
following conditions: 

1. ^(0) = 1 

2. ipioo) = limt^c>oip{t) = 

3. ijj is strictly decreasing on [0, inf{t : ip{t) = 0}] 
Definition: Archimedean Copula 

A d- dimensional copula C is called Archimedean if for some generator i/j it can be represented as: 

C(u) = ^{^-i(ni) + --- + ^-Hnrf)} =V{i(u)} Vu G [0, l]'^ (4.2) 

where : [0, 1] — t- [0,oo] is the inverse generator with ip~^{0) = inf{t : ip{t) = 0}. 

Note the shorthand notation t(u) = ^~^(ui) + • • • + '(p~^{U(i) that will be used throughout this section. 

As we will see later, it is necessary to have formulas for computing the copula densities if one seeks 
to fit these models using a maximum likelihood approach. Equation 14.31 provides such a formula in a 
generic form for each member of the family of Archimedean copulae. 

Definition: Archimedean Copula Density 



McNeil and Neslehovd 12009 ) prove that an Archimedean copula C admits a density c if and only ifip^^^^) 
exists and is absolutely continuous on (0,oo). When this condition is satisfied, the copula density c is 
given by 

c(u) = ^'27.\'.du'''^ = V'^'^i^HinCV'-'/C^.) > ^^(0,1)'^ (4.3) 
1 ■ ■ ■ j=l 



4.1.1 One-parameter Archimedean Members: 

In this section we describe three of the one parameter multivariate Archimedean family copula models 
which have become popular model choices and are widely used for estimation. This is primarily due to 
there directly interpretable features. We select these three component members, the Clayton, Frank and 
Gumbel models, for our mixture models since they each contain differing tail dependence characteristics. 
Clayton provides lower tail dependence whilst Gumbel provides upper tail dependence. The Frank copula 
also provides dependence in the unit cube with elliptical contours with semi-major axis oriented at either 
7r/4 or 37r/4 depending on the sign of the copula parameter in the estimation. Therefore the Frank model 
component will allow us to capture parsimoniously potential negative dependence relationships between 
the currencies in the portfolio under study. Formulas for these copulae, as well as their respective gener- 
ators, inverse generators and the d-th derivatives of their generators (required for the density evaluation) 
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are given i n Table [H The exp licit formulas for the d-th derivatives for all of the copulae in Tabled] were 
derived in Hofert et alJ pOld ). 



4.1.2 Two-parameter Archimedean Members via Outer Power Transforms 



In this section we also consider more flexible generalizations of the single parameter Archimedean members 
discussed above. To achieve these generalizations we consider the outer-power transforms o f the Clayton, 
Prank and Gumble members, as discussed below which is based on a result in iFelleii (jl97ll ). 
Definition: Outer power copula 

i_ 

The copula family generated by \lj{t) = ip{ti^) is called an outer power family, where /3 G [l,oo) and 
V' £ ^oo (the class of completely rnonotone Ar chimedean generators). 



The proof of this follows from 



Felled (|l97lh . i.e. the composition of a completely monotone function 



with a non-negative function that has a completely monotone de rivative is again completely monotone. 
Such copula model transforms were also studied in iNelsenI (jl997l ). where they are referred to as a beta 



family associated with the inverse generator ip~^. 

As has been noted above, in performing the estimation of these transformed copula models via likeli- 
hood based inference it will be of great benefit to be capable of performing evaluation pointwise of the 
copula densities. In the case of the outer power transformed models, this will require the utilization of 
a spec ific multivariate chain rule resu 



(j 18571 ) and discussions in for example 



t widely known as the Faa di Bru no's Forniula, see 



Faa di Bruno 



Constantine and Savitd ()1996l ) and iRomanI (j 19801 ) . To understand 



how such a result is required consider the following remark. 

Remark: The generator derivatives for the outer power transforms can be calculated using the base 
generator derivatives and the following multi-dimensional extension to the chain rule for the outer power 
versions. The densities for the outer power copulae in TableU\can thus be calculated using equation \4.3[ 
Before stating Faa di Bruno's Formula for differentiation of multivariate composite functions via a 
generalized chain rule, it will be convenient notationally to present such results with respect to Bell 
polynomials. Therefore w e recall t h e def inition of such polynomials below, which are widely used in 
combinatorics analysis, see 



Mihoubi 



(20081) for details. 



Definition: Bell Polynomial 

The Bell polynomial with arguments n and k is given by 



Bn,k{xi,X2, ... ,x. 



n-fe+l ) 



E 



ni 



Xn-k+l 



(4.4) 

■ = k 



{n-k + l)\ 

where the sum is taken over all sequences Ji, J2> Jn-fc+i of non-negative integers such that ji +j2 + 

and ji + 2j2 + 3j3 H = n. 

These polynomials are then utilised to simplify the expressions for the differentiation of multivariate 
composite functions in Faa di Bruno's Formula as detailed next. 
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(4.5) 



Faa di Bruno's Formula: iRiordanI (j 19461 ) 

If f and g are functions with a sufficient number of derivatives, then 

k=0 

where B^^k o-'''^ the Bell polynomials, defined below. 

Having specified the details of the parametric models from the Archimedean family that we are con- 
sidering, it is also instructive to understand from the perspective of their dependence properties why we 
are interested in considering such models. To address this question, it is typically the case that one will 
resort to a discussion on what is known as tail-dependence. 



4.2 Tail Dependence 

In order to examine the dependence behaviour in the extremes of the multivariate distribution we use 
the concept of tail dependence. The bivariate tail dependence coefficient is defined as the conditional 
probability that a random variable exceeds a certain threshold given that the other random variable in the 
jution has exceeded thi s thres hold. This concept has been recently extended to the multivariate 



joint distri 

setting bv lDe Luca and Riviecciol (|2012l ). Now one may accurately interpret the tail dependence present 
between sub- vector partitions of the multivariate random vector with regard to joint tail dependence 
behaviours. In the context of the applications we consider in this paper, this allows us to examine the 
probability that any subvector of the log return forward exchange rates for the basket of currencies will 
exceed a certain threshold given that the log return forward exchange rates for the remaining currencies 
in the basket have exceeded this threshold, in particular thresholds that are placing an interest in the 
tails of the multivariate distribution. The interpretation of such results is then directly relevant to 
assessing the chance of large adverse movements in multiple currencies which could potentially increase 
the risk associated with currency carry trade strategies significantly, compared to risk measures which 
only consider the marginal behaviour in each individual currency. 

Below we give the definitions of upper and lower tail dependence, as well as the e xplicit generalised mul 



tivariate expressions for Archimedean copulae, equations 14.61 and 14.71 derived in 



De Luca and Rivieccio 



Definition: Generalized Upper Tail Dependence 

Let X = {Xi, ...,Xd)^ be a d dimensional random vector with marginal distribution functions Fi, ...,Fd. 
The coefficient of upper tail dependence is defined as: 



XyMh+i,...,d ^ lim P{X, > F-Hu),...,Xh > F-Hu)\Xh+i > F~\v),...,Xd > F-\v)) 



lim 



Eti( 




i)~^'{it) 


) 


Tt^ (( 





(4.6) 
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Definition: Generalized Lower Tail Dependence 

Let X = {Xi, ...,X^)'^ be a d dimensional random vector with marginal distribution functions Fi, ...jF^. 
The coefficient of lower tail dependence is defined as: 

^i,...Mh+i,...,d ^ lim P{X, < F-Hu),...,Xh < < X, < F-^u)) 



d ij^^'idt) 



(4-7) 



lim 

t^ood-h'llj-^'{{d-h)t) 



4.2.1 Mixtures of Archimedean Copulae 



In order to add additional flexibility in the possible dependence features one can study for the currency 
portfolios, we decided to utilize mixtures of copula models. In this regard we have the advantage that we 
can consider asymmetric dependence relationships in the upper tails and the lower tails in the multivariate 
model. In addition we can perform a type of model selection purely by incorporating into the estimation 
the mixture weights associated with each dependence hypothesis. That is the data can be utilised to 
decide the strength of each dependence feature as interpreted directly through the estimated mixture 
weight attributed to the feature encoded in the particular mixture component from the Archimedean 
family. 

In particular we have noted that mixture copulae can be used to model asymmetric tail dependence, i.e. 
by combining the one-parameter or two-parameter families discussed above or indeed by any combination 
of copulae. This is possible sin ce a linear convex combination of 2 copulae is itself a copula, see discussions 



on this result in 



NelsenI pOOd ). 



Definition: Mixture Copula 

A mixture copula is a linear weighted combination of copulae of the form: 



N 



i=l 



where 0<Aj<l \/i = 1, N and X^i^i — ^ 

Thus we can combine a copula with lower tail dependence, a copula with positive and negative depen- 
dence and a copula with upper tail dependence to produce a more flexible copula capable of modelling 
the multivariate log returns of forward exchange rates of a basket of currencies. For this reason in this 
analysis we will use the Clayton-Frank-Gumbel mixture model. In addition to the C-F-G mixture model 
we will also investigate a mixture of outer power versions of the base copula Clayton, Frank and Gumbel. 
Remark We note that the tail dependence of a mixture copula can be obtained as the linear weighted 
combination of the tail dependence of each c ompo nent i n the mixture weigh ted by the appropriate mixture 



weight, as discussed in for example 



NelseA IMM) and 



Peters et al 



(20h 
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Table 1: Archimedean copula generator functions, inverse generator functions and generator function d-th derivatives. 



Family 



Clayton 



OP-Clayton 



1 +t!i 



{s-P - 1) 
{s-P - If 



13" 



Frank 



OP-Frank 



iln [l - e-*(l - e 



l-e-'^{l-e 



In 



In 



itLi_(,_i){(l - e-P)e-'] 



1 \ fe 



Gumbel 



OP-Gumbel 



1 

-tp 



(-Ins) 



(-Ins) 



Rem ark: The densities for the one-parameter copulae in Tabled] can be calculated using equation 14.31 For details of the results contained in this 



table see 



Hofert et al 



(12013) 



= HE ■=!(■) - k e i,...,d 
^ Lis{z) = j2T=i w 



5 Likelihood Based Estimation of the Mixture Copula Models 

We begin this section with a discussion on the choices we make for the marginal distributions for each 
of the currencies specified in the baskets constructed for the high interest rate differentials and also the 
baskets for the low interest rate differentials. 

In modelling parametrically the marginal features of the log return forward exchange rates, we wanted 
flexibility to capture a broad range of skew-kurtosis relationships as well as potential for sub-exponential 
heavy tailed features. In addition, we wished to keep the models to a selection which is efficient to 
perform inference and easily interpretable. We therefore considered a first analysis utilizing LogNormal 
distributions for the monthly forward exchange rate returns, which would be equivalent to specification 
of a Normality assumption on the distribution for the log return forward exchange rates. This model is 
given by the following parameteric density, for a random variable X ~ F[x] /i, o"), in Equation 15.11 below. 

with the shape parameter cr^ > and the log-scale parameter /i G M and the support x G (0, oo). 

We found when analysing the goodness-of-fit for this LogNormal model on each of the assets in the 
20 currencies considered, over both 6 month and 1 year sliding windows, that the fit of the LogNormal 
model would be systematically rejected as a suitable model for a couple of currencies. In the majority 
of cases over these sliding windows (locally stationary time series) the LogNormal model was more than 
adequate. However, since some of the currencies that were rejecting this fit were appearing regularly in 



the high interest rate baskets we also decided to consider a more flexible three parameter mode^ 
margi nal d istributions given by the Log-Generalized-Gamma distribution (l.g.g.d.), see details in 



for the 



Lawless 



(119801) and lGonsul and Jainl (jl97l|). 

The l.g.g.d. is a parametric model based on the generalized gamma distribution which is highly utilized 
in lifetime modelling and survival analysis. The density for the generalized gamma distribution and the 



l.g.g.d are given respectively by Equations 15.21 and 15.31 

/x(x;^,a,/3) = ^^exp^-y J (5.2) 

with parameter ranges /c > 0, a > and /3 > and a support of x G (0, oo). Then the log transformed 
g.g.d. random variable 1" = InX is given by the density of the l.g.g.d. as follows. 



fYiy;k,u,b) = ^p^exp 



k — - — — exp 



(5.3) 



h J \ b 

with u = log(a), b = (]~^ and the support of the l.g.g.d. distribution is y G M. 

This more flexible three parameter model is particularly interesting in the context of the marginal 
modelling we are considering since the LogNormal model is nested within the g.g.d. family as a limiting 
case. In addition the g.g.d. also includes the exponential model {(3 = k = 1), the Weibul distribution with 
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(k = 1) and the Gamma distribution with (/3 = 1). Next we discuss how one can perform inference for 
the multivariate currency basket models using these marginal models and the mixture copula discussed 
previously. 



5.1 Two Stages: Inference For the Margins 



The inference function for margins (IFM) technique introduced in lJod (j2005l ) provides a computationally 
faster method for estimating parameters than Full Maximum Likelihood, i.e. simultaneously maximising 
all model parameters and produces in many cases a more stable likelihood estimation procedure. An 
alternative approach to copula model parameter estimation that is popular in the lit erature is known as 



Genest et al. 



the Maximum Partial Likelihood Estimator (MPLE) detailed in 

The procedure we adopt for likelihood based estimation is the two stage estimation known as Inference 
on the Margins which is studied with regard to the asymptotic rel ative efficie ncy of t he two-stage estima 



Joel (120051) and in 



Hafner and Manner 



tion p rocedure compared with maximum likelihood estimation in 
(|2010l ). It can be shown that the IFM estimator is consistent under weak regularity conditions. How- 
ever, it is not fully efficient for the copula parameters. Nevertheless, it is widely used for its ease of 
implementation and efficiency in large data settings such as the models we consider in this study. 

To complete this discussion on general IFM, before providing the MLE estimation expressions, we first 
note that in our study we fit copula models to the high interest rate (IR) basket and the low IR basket 
updated for each day in the period 04/01/2000 to 02/01/2013 using log return forward exchange rates at 
one month maturities for data covering both the previous 6 months and previous year as a sliding window 
analysis on each trading day in this period. Next we discuss briefiy the marginal MLE estimations for 
the LogNormal and the l.g.g.d. models. 



5.1.1 Stage 1: Fitting the Marginal Distributions via MLE 

In the first step we fit the marginal distributions to either the LogNormal model or the l.g.g.d model. In 
the case of the LogNormal model this is achieved effortlessly since we may utilise the well-known analytic 
expressions for the MLE estimates: 

^^log(x,) 

(5.4) 



In the case of the l.g.g.d. distribution the estimation for the three model parameters can be significantly 
more challenging due to the fact that a wide ran ge of model p arameters, especially for k can produce 



similar resulting density shapes, see discussions in 



Lawless! (jl98d ). To overcome this complication and to 



make the estimation efficient it is proposed to utilise a combination of profile likelihood methods over a 
grid of values for k and perform profile likelihood based MLE estimation for each value of k, then for the 



19 



other two parameters b and u. The differentiation of the profile likehhood for a given value of k produces 
the system of two equations given by 

cr\/fc 

1 ^ / „, \ 

exp(/i) 



1 " / 



1=1 

Vi \ 



(5.5) 



y 



with n the number of observations, i/i = logXj and the parameter transformations ^ = ^^'^ = 
u + hlnk. The second equation is solved directly via a simple root search for the estimation of a and 
then substitution into the first equation provides the estimation of jl. Note, for each value of k we select 
in the grid, we get the pair of parameter estimates fx and ct, which can then be plugged back into the 
profile likelihood to make it purely a function of /c, with the estimator for k then selected as the one with 
the maximum likelihood score. 

5.1.2 Stage 2: Fitting the Mixture Copula via MLE 

In order to fit the Clayton-Frank-Gumbel model the copulae parameters {pdaytom P frank, Pgumbei) and 
the copulae mixture parameters {Xdayton, ^frank, ^gumbei) are estimated using maximum likelihood on the 
data after conditioning on the selected marginal distribution models and their corresponding estimated 
parameters obtained in Stage 1. These models are utilised to transform the data using the cdf function 
with the mle parameters {jl and a) if the LogNormal model is used or (A;, u and h) if the l.g.g.d is 
considered. 

Therefore, in this second stage of MLE estimation we aim to estimate either the one parameter 
mixture of C-F-G components with parameters e_ = {pdayton, P frank, PgumbeuKiayton A frank, \umbei) or 
the two parameter mixture of outer power transformed mixture components OC-OF-OG components 

with parameters 0_ — {Pdayton, Pfrank, Pgumbel, ^dayton, ^frank, Xgumbel, Pdaytom Pfrank, l^gumbeO- This is 

achieved in each case by the conditional maximum likelihood. To achieve this we need to maximise the 
log likelihood expressions for the mixture copula models, which in our framework are given generically 
by following function for which we need to find the mode, 

n n d 

m = J2^''9 c^~^'^{Fi{Xa;pi,ai),...,Fd{Xid;pd,^d)) + ^^^log fj{X,j; p,,aj) (5.6) 

i=l 1=1 j=l 

with respect to the parameter vector 9. 

For example in the case of the Clayton-Frank-Gumbel mixture copula we need to maximise on the 
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log-scale the following expression. 

n 

= [Ac * {cpc {Fi (Xa; Ai, <Ti) ...,Fd (X^^; /i^, a^))) 

i=l 

(5.7) 

+ Xf* {c^p [Fi {Xii;fii,ai) . . . , Fd{Xid; fid,^d))) 

+ ^G* {c^a iFiiXii;fii,ai) ...,Fd{Xid;fid,^d)))] 

This optimization is achieved via a gradient descent iterative algorithm which was found to be quite 
robust given the likelihood surfaces considered in these models with the real data. To illustrate this 
point, at this stage it is instructive to present some examples of the shapes of the profile likelihoods 
that are being optimized over for some of the important copula model parameters in the C-F-G mixture 
example for a 6 month window of data randomly selected from the data set for both the high interest 
rate basket and the low interest rate basket. Example plots of the profile likelihood for the 6-dimensional 
optimisation space for two different example days can be seen in Figures [3] and [H 

5.2 Goodness-of-Fit Tests 

In this section we briefly comment on the model selection aspects of the analysis we undertook. As 
mentioned we flrst undertook a process of fitting the marginal LogNormal model to all of the 20 currencies 
considered in the analysis over a sliding window of 6 months and 1 year. For each of these fits we then 
performed a formal hypothesis test in which we postulated that the null distribution is the LogNormal 
model and then look for evidence in the data to reject this hypothesis at a level of significance of 5%. 
To undertake this test we considered the standard Kolmogorov-Smirnov test. As we will present in 
the results we found strong evidence to reject the null systematically for a few important developing 
countries' marginal models, hence we also undertook estimation of the l.g.g.d. models for all of the 20 
currencies. We are particularly interested in this case in the optimal choice of the model parameter k 
which as it asymptotically gets large A; — )■ oo will produce a LogNormal model. We found as expected the 
estimated model fits were significantly improved when fitting the l.g.g.d. models for the cases in which 
the LogNormal was rejected by the K-S test. In addition the estimated k parameter in the periods of 
rejection of the LogNormal hypothesis were estimated at values significantly lower than the upper bound 
in our search space. We assessed the optimal choice of marginal model between the LogNormal and 
the l.g.g.d. models then via a standard information criterion based on the Akaike Information Criterion 
(AIC). 

In terms of the selection of the copula mixture models, between the mixture of one parameter C-F-G 
model versus the two parameter mixtures of OC-OF-OG models, we again used a scoring via the AIC. 
We note that there are also alternative information criterion developed for copula models to assess the 
joint suitability of the copula model incorporating both the marginal and the joint copula structure which 
are modifications of the AIC, adjusting the penalty term for the approach adopted in the estimation, see 



21 



for example the Copula-Information-Criterion (CIC) in lGr0nneberd ((20101) for details. The results are 
presented for this comparison in Figure [5] in the top panel for the high interest rate basket and in Figure 
[5] in the lower panel for the low interest rate basket, over time based on the 6 month sliding window. 

To further analyse this comparison of optimal copula mixtures we plot the AIC differentials for each 
of the currency baskets in Figure [6l 

Figures and E] show it is not unreasonable to consider the C-F-G model for this analysis, since the 
mean difference between the two AIC scores for the models is 2.05 in favour of the C-F-G. However, we 



do note that the OP.C-OP.F-G model seems to fit better during crisis periods. 



6 Results 

In this section we present a detailed analysis of the estimation of the marginal distributional models 
and the mixture copula models for both the high interest rate basket and the low interest rate basket. 
Firstly, we investigate the properties of the marginal distributions of the exchange rate log-returns for 
the 20 currencies. We then interpret the time- varying dependence characteristics of the fitted copula 
models to the high interest rate basket and the low interest rate basket across the period 04/01/2000 to 
02/01/2013. Note, all results presented below are for the case in which we considered a 6 month sliding 
window, results for the 1 year sliding window were similar in nature and so are omitted. 

6.1 Modelling the Marginal Exchange Rate Log- Returns 

In order to model the marginal exchange rate log-returns we first fit LogNormal models to each of the 20 
currencies considered in the analysis, updating the fits for every trading day in the period 04/01/2000 to 
02/01/2013 based on the 6 months sliding window. The LogNormal model was selected due to the fact it 
has a positive support, represents a range of skew-kurtosis characteristics and can display sub-exponential 
tail features (ie. heavy tailed features) should such attributes be present in the data. We assessed the 
quality of the fits for each currency using a standard Kolmogorov-Smirnov goodness-of-fit test, at the 
5% significance level. A summary of the results of this analysis are presented in Table [2] which shows the 
proportion of rejections of the null hypothesis, that the marginal distribution is LogNormal for each of 
the currencies on a given 6 months block of trading days. 

We learn from this analysis that the majority of the currencies demonstrate reasonable marginal 
distribution fits under a LogNormal family, however there are a few notable exceptions. Specifically the 
Turkish Lira, Malaysian Ringgit, Indian Rupee, Thai Baht, South Korean Won and Taiwanese Dollar 
demonstrated sustained periods in the analysis in which the LogNormal model would be unsuitable to 
capture the features of the time series adequately. This is significant in this analysis since these currencies 
actually correspond to the currencies that have a strong presence in the high interest rate baskets, as seen 
in Figure [TJ Therefore, they will play an important role in the multivariate analysis of the currency carry 
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trade. As such, it is important to accurately model the features of each of these particular currencies' 
marginal distributions, before undertaking the multivariate mixture copula analysis, we proposed to 
generalize the marginal model analysis to a more flexible three parameter family of models given by the 
log generalized gamma distribution, as discussed in Section 5. 

The log-generalised gamma distribution (l.g.g.d.) should improve the fit for all currencies since it allows 
for more flexibility in the tails of the distribution and a wider range of skew-kurtosis relationships when 
compared to the LogNormal model family. In addition, as we note in Section 5, for those currencies in 
which the LogNormal model was a suitable fit, then they will still obtain such distributional characteristics 
since the LogNormal model is a limiting case of the l.g.g.d. as k tends to infinity. Hence, we can still 
incorporate the LogNormal model for the currencies that were a good fit. 

The maximum likelihood parameters {fi,a,k) of the fitted l.g.g.d. margins for each of the currencies 
can be seen in the Appendix in Figures [TJ [5] and O These plots demonstrate the time varying attributes 
of the marginal distributions for each currency, illustrating interesting changes in tail behaviour and 
skewness-kurtosis characteristics over time, especially in heightened periods of volatility in some of these 
currencies. In particular, there are three standout periods (2003, 2009 and 2012) of heightened fi and 
a parameter values across most of the currencies. Hence, during these periods the exchange rate log- 
returns may demonstrate heavier tails, and increased volatility in the parameter estimates. In addition, 
we observe that a few important currencies for the currency carry trade analysis demonstrate sustained 
differences in their marginal distribution attributes relative to the other currencies. An important example 
of this is the fi estimates in Figure [7] for the TRY, the NZD and the BRL. Similar significant differences 
between these particular currencies and the rest of the currencies are observed in the estimates of a in 
Figure [HI 

As the value of the parameter k in the l.g.g.d. gets large we expect the LogNormal fit to be a suitable 
model structure for the marginal distributions. As illustrated in the K-S test results certain currencies 
systematically did not have a suitable fit with the LogNormal model. Examples of this are clear when 
we consider the estimates of k in Figure O Again we see systematically smaller values for the estimate 
of k in the TRY and the BRL. We see clearly in Figure [9] the periods of time during which the currencies 
display non LogNormal behaviour. The most prominent example being the Turkish Lira (orange), which 
shows consistently low values of A;. As we noted in Section 5, for small values of A; ~ 1 we obtain WeibuU 
like tail behaviour and in addition, in the cases when a ~ 1 jointly with small values of k, we expect the 
light tailed exponential models to be suitable. As a consequence of this analysis and comparison of AIC 
results we proceeded with the joint estimation utilising the l.g.g.d. marginal models for every currency. 

A noticeable period for the Turkish Lira is early in 2001 during which low values of the parameter 
k clearly provides evidence of heavy tail log-returns distribution for this specific currency. We indeed 
mentioned earlier in this paper the Turkish government's decision in February 2001 to stop draining 
reserves to bolster its currency led the same day to a 30% devaluation of the Turkish Lira relative to the 
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dollar. 

Table 2: Proportion of rejections of the null hypothesis that the sample is from a LogNormal 
distribution, measured using a k-s test at the 5% level. 

Block length EUR TRY JPY GBP AUD CAD NOK CHF SEK MXN 

6 month 0.001 0.198 0.043 0.000 0.023 0.000 0.000 0.031 0.012 0.032 
Year 0.000 0.553 0.107 0.007 0.120 0.018 0.006 0.084 0.018 0.128 



Block length PLN MYR SGD INR ZAR NZD THB KRW TWD BRL 

6 month 0.018 0.494 0.000 0.234 0.025 0.012 0.221 0.130 0.192 0.086 
Year 0.094 0.651 0.071 0.549 0.124 0.113 0.504 0.350 0.381 0.403 

6.2 Copula Modelling Results 

We now utilised each of the l.g.g.d. marginal distribution fits for a given day's set of currencies in the 
high interest rate and low interest rate baskets to analyse the joint multivariate features. To achieve this 
for each of the currencies, the exchange rate log-return data was transformed via the l.g.g.d. marginal 
model's distribution function to uniform [0, 1] margins. Then the mixture Clayton-Frank-Gumbel copula 
(denoted C-F-G ) and the outer-power versions were fitted each day to a sliding window of 6 months and 
one year log-returns data for both the high interest rate and low interest rate baskets. Below we will 
examine the time- varying parameters of the maximum likelihood fits of this mixture C-F-G copula model. 
Furthermore, the results for the outer-power transform cases did not demonstrate discernible differences 
from the base C-F-G model and so were excluded. This can be seen from the figures displaying the AIC 
for each of these models. 

In this analysis there are several attributes to be considered for the mixture copula model, such as the 
relevant copula structures for the high and low interest rate baskets and how these copula dependence 
structures may change over time. In addition, there is the strength of the tail dependence in each currency 
basket and how this changes over time, especially in periods of heightened market volatility. The first 
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of these attributes we will consider to be a structure analysis studying the relevant forms of dependence 
in the currency baskets and the second of these attributes that we shall study will be the strength of 
dependence present in the currency baskets, given the particular copula structures in the mixture. 

Therefore we first consider the structural components of the multivariate copula model. To achieve 
this, we begin with a form of model selection in a mixture context, in which we consider the estimated 
relative contributions of each of the copula components (and their associated dependence features) to 
the joint relationship in the high and low interest rate currency baskets over time. This is reflected in 
the estimated mixture component weights, which can be seen in Figures [10] and [11] for the high interest 
rate basket and low interest rate basket respectively. The A values show the relevance of each of the 
component copulae to the data. Thus a small A value indicates the lack of a need for that particular 
copula component in order to model the associated 6 months or one year block of data. In contrast, 
for example a A value for the Gumbel component very close to 1 indicates the block of data could be 
well modelled by a Gumbel copula alone. Hence, these plots convey the time varying significance of 
hypotheses about the presence of upper and lower tail dependence in each of the baskets over time. 
Examining these plots shows that in general the Clayton mixture weight tends to be lower when the 
Gumbel mixture weight is higher. We can also see that the Frank copula is systematically present in 
the mixture. In addition, we see that in the periods of high market volatility we observed differences 
in the relevant upper and lower tail dependence structural attributes when comparing the high versus 
low interest rate baskets. That is, there is an asymmetric tendency for the presence of particular copula 
components over time when comparing the high and low interest rate baskets. The implications of this 
will be discussed in further detail in the discussions. 

In terms of the second attribute, the strength of the copula dependence, we analyse this in several 
ways. Firstly through an analysis of the estimation copula parameter components over time, then through 
an analysis of the transformation of these copula parameters to rank correlations and finally through an 
analysis of the multivariate strength of the mixture copula tail dependence over time. 

The individual component copula parameters can be seen in Figures [T2] and [13] for the high interest 
rate basket and low interest rate basket respectively. The strength of the copula parameters in the 
baskets shows a large degree of variance during the period 04/01/2000 to 02/01/2013. One interesting 
observation is the very large spikes in the Gumbel copula parameter observed for the high interest rate 
basket from 2006 to 2007 and again in 2009. This was significant as it also corresponds to periods in 
which the Gumbel copula mixture weight was non-trivial. 

The measure of concordance as captured by Kendall's tau is decomposed in this analysis according 
to each of the mixture components, scaled by the mixture weights A, and can be seen in Figure [TJ] for 
the high interest rate basket and Figure [TCI for the low interest rate basket. These plots provide a more 
intuitive picture of the time-varying contributions of the individual copulae to the dependence structure 
present in each of the baskets. Interestingly, we see the rank correlation contribution from the Frank 
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copula indicates the presence of negative as well as positive rank correlations. In addition, as discussed 
with the mixture weights, there is perhaps some asymmetry present between the high and low interest 
rate baskets over time. 

Perhaps the most interesting and revealing representation of the tail dependence characteristics of 
the currency baskets can be seen in Figures [16] - EH Here we can see that there are indeed periods 
of heightened upper and lower tail dependence in the high interest rate basket. There is a noticeable 
increase in upper tail dependence at times of global FX volatility. Specifically, during late 2007, i.e. the 
global financial crisis, there is a sharp peak in upper tail dependence. Preceding this, there is an extended 
period of heightened lower tail dependence from 2004 to 2007, which could tie in with the building of the 
leveraged carry trade portfolio positions. 

In understanding this analysis we note that Figures [16] and [T7] show the probability that one currency 
in the basket will have a move above/below a certain extreme threshold given that the other three 
currencies have had a move beyond this threshold. Then in Figures [18] and [19] we show the probability 
that two currencies in the basket will have a move above/below such an extreme threshold given that the 
other two currencies have had a move beyond this threshold. Finally, in Figures [20] and \TI] we show the 
probability that three currencies in the basket will have a move above/below a certain threshold given 
that the remaining currency has had a move beyond this threshold. 

To illustrate the relationship between heightened periods of significant upper and lower tail dependence 
features over time and to motivate the clear asymmetry present in the upper and lower tail dependence 
features between the high and low interest rate baskets over time we consider a further analysis. In 
particular we compare in Figures [22] and [23] the tail dependence plotted against the VIX volatility 
index for the high interest rate basket and the low interest rate basket respectively for the period under 
investigation. The VIX is a popular measure of the implied volatility of S&P 500 index options - often 
referred to as the fear index. As such it is one measure of the market's expectations of stock market 
volatility over the next 30 days. We can clearly see here that in the high interest rate basket there are 
upper tail dependence peaks at times when there is increased stock market volatility, particularly post- 
crisis. However, we would not expect the two to match exactly since the VIX is not a direct measure of 
global FX volatility. We can thus conclude that investors' risk aversion clearly plays an important role in 
the tail behaviour of high interest rate currencies and more importantly in their dependence structure. 
This statement can also be associated to the globalization of financial markets and the resulting increase 
of the contagion risk between countries. This conclusion corroborates some of the recent literature results 
with regards to the skewness and the kurtosis features characterizing the c urrency carry trade portfolios 
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The black lines plotted in Figures [25] and [23 furthermore display the mean tail dependence before and 
after August 2007 (which corresponds to the beginning of the global financial crisis). The data shows 
a large increase in upper tail dependence in the high interest rate basket after the crisis, as well as a 
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smaller decrease in lower tail dependence. Interestingly there is very little difference in the mean tail 
dependence before and after the crisis for the low interest rate basket. The carry trade portfolios were 
particularly impacted by the sub-prime crisis as most of these currency positions were implemented and 
held by financial institutions which faced sudden difficulties to finance the leverage of their positions. 
Furthermore, another interesting point we can make from the analysis of these two figures is the higher 
level of lower tail dependence before the financing crisis, especially between 2004 and 2007. If we put in 
parallel the fact that during this three year period the VIX index was noticeably and continuously going 
down we could imagine that this increase of the lower tail dependence results from lower risk aversion 
and the resulting tendency of investors to accordingly increase their leverage on risky positions such as 
currency carry trades. 

7 Discussion 

This article sheds light on an important topic that has been discussed in the recent financial literature 
which is the forward premium puzzle. This market phenomenon is particularly interesting from a theoret- 
ical standpoint as well as for the understanding of financial market mechanisms. It has been demonstrated 
that the currency markets were indeed not respecting empirically a fundamental relation in finance which 
connects the currency exchange rates and the interest rates associated with two different countries. The 
main contribution of this article has been to propose a rigorous statistical modelling approach which 
captures the specific statistical features of both the individual currency log-return distributions as well 
as the joint features such as the dependence structures prevailing between all the exchange rates. 

In achieving this goal, we first assessed the marginal statistical features of each of the 20 currencies 
on an assumed locally stationary sliding window of six months, over all the trading days in the period 
04/01/2000 to 02/01/2013. We found that a simple LogNormal marginal distribution would not produce 
a suitable statistical fit for some of the key currencies that are regularly present in the high interest rate 
basket throughout this period. As detailed in the results section this was notably the case in unstable 
economies such as developing countries (for instance Turkey, Brazil or South Africa) where political 
stability or default risk create sudden and violent adjustments to their currency exchange rates with 
other countries. We note that these currencies are still of direct significance to currency carry trade 
strategies since according to the modern portfolio theory this intrinsic risk borne by an investor in these 
currencies can be diversified and mitigated by adding to the considered portfolio other currencies which 
depend themselves on different sources of intrinsic risk. This would effectively establish a diversified 
portfolio of currencies violating the UIP hypothesis and would thus provide a very attractive average 
return for a very limited risk which has been the conclusion of several recent empirical studies in the 
finance literature. 

The conclusion of this is that we cannot exclude these currencies from our high interest rate basket 
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analysis, even though they may demonstrate attributes resulting primarily from significant changes in 
their countries political and financial structure. As a result we needed to obtain more flexible marginal 
models to capture the features of these currencies more adequately. Consequently we modelled each 
currency exchange rate return marginally via a flexible three parameter parametric model which offers a 
wide range of skew-kurtosis relationships as well as the possibility of light exponential tails and heavier 
sub-exponential tail behaviours such as the LogNormal member. The parametric family of distributions 
we selected for this purpose was the log-generalized gamma distribution. 

Having modelled the marginal attributes of the high and low interest rate currency baskets over 
time adequately, it was then our main emphasis to assess the multivariate dependence features of the 
currency baskets. In particular how this may change over time within a given basket, where we were 
particularly interested in the effect of the composition of the basket over time, and the response of the 
multivariate dependence features of the modelled basket and how it may respond in periods of heightened 
market volatility versus more stable periods. In addition to this within basket temporal analysis, from 
the perspective of undertaking a currency carry trade strategy, we would need to consider the relative 
relationships between the temporal dependence features of the high interest rate and low interest rate 
currency baskets. We demonstrate several interesting features from our model fits relating to asymmetries 
between the high and low interest rate baskets over time, especially during periods of high volatility in 
global markets. One way we ascertained such periods was through a comparison of the VIX versus features 
of the multivariate dependence relationships we modelled. Importantly we found substantial evidence to 
support arguments for time varying behaviours both in the structural dependence hypotheses posed about 
the currency baskets, as captured by the relevant contributing copula components to the multivariate 
mixture model. As well as substantial evidence for significant tail dependence features in both the high 
and low interest rate baskets, which again displayed interesting asymmetries between the high and low 
interest rate baskets over time. 

The financial interpretation of the significance of these findings is related to the fact that it demon- 
strates that historically average rewards from a currency carry trade portfolio can be exposed to a 
significant risk of large losses arising from joint adverse movements in the currencies that would typically 
comprise the high and low interest rate baskets that an investor would go long and short on when trading. 
Hence, we conclude that our second contribution to the literature has been to rigorously demonstrate 
that such assertions relating to the profitability of the currency carry trade are failing to appropriately 
take into consideration an important component of the risk which characterizes these types of portfolios 
of currencies named carry trade portfolios. 

We conclude that indeed the copula theory we employed allows us to demonstrate statistically that 
beyond the intrinsic risk associated to high interest rate countries (which are generally paying higher 
interest rates to compensate for a higher risk) typically studied in the literature from a marginal per- 
spective, another source of risk plays an important role. This second source of risk is related to the 
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dependence structures linking these high interest rate currencies, more specifically the significant tail de- 
pendence features observed in our model analysis. We indeed proved through a mixture of Archimedean 
copulae the significant presence of tail dependence among high interest rate currencies which could have 
dramatic consequences on the carry trade portfolio's risk profile when accounted for appropriately in risk 
reward analysis. As a matter of fact, the tail dependence directly influences the diversity of the assets 
and thus reduces the appealing convergence property stated by the modern portfolio theory. 

Said differently, our copula based probabilistic modelling approach allows us to demonstrate that be- 
sides the intrinsic risk associated to each particular high interest rate currency, another factor constitutes 
a determining source of risk which turns out to be the level of risk aversion prevailing in the market. It 
was demonstrated in our analysis that both upper and lower tail dependence features displayed signifi- 
cant association and asymmetries with each other between the high and low interest rate baskets during 
periods of relative financial stability versus periods of heightened market volatility. 

These tail dependence features in the high interest rate basket were significantly increasing during crisis 
periods leading to an increased amount of risk associated with utilising such currency baskets (which were 
no longer diversified due to the presence of significant tail dependence features) in a carry trade. That 
being said, a rational portfolio manager's natural risk aversion tells them that they should receive an 
additional remuneration in order to offset any additional sources of risk associated to an investment. 
Therefore, to properly assess the profitability of the currency carry trade, such tail dependence features 
should be incorporated into the analysis of such risk-rewards when developing a trading strategy. To 
conclude, our paper rigorously tempers the too often claimed attractiveness of the currency carry trade 
and provides to investors a risk management tool in order to control and monitor the risk contained in 
such positions. 
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Figure 3: Profile likelihood plots for C-F-G mixture model. 
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Figure 4: Profile likelihood plots for C-F-G mixture model. 
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Figure 5: AIC comparison of C-F-G vs OP.C-OP.F-G for 6 month blocks on high and low 
IR baskets. 
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Figure 6: AIC comparison of C-F-G vs OP.C-OP.F-G for 6 month blocks on high and low 
IR baskets. 
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Figure 7: ^ parameter of log generalised gamma margins using 6 month blocks 
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Figure 9: k parameter of log generalised gamma margins using 6 month blocks 
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Figure 10: A Mixing proportions of the respective clayton, frank and gumbel copulas on the 
high interest rate basket, using 6 month blocks. 
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Figure 11: A Mixing proportions of the respective clayton, frank and gumbel copulas on the 
low interest rate basket, using 6 month blocks. 
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Copula parameters using 6 month blocks for high IR basket 
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Figure 12: p Copula parameters for the clayton, frank and gumbel copulas on the high interest 
rate basket, using 6 month blocks. 
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Figure 13: p Copula parameters for the clayton, frank and gumbel copulas on the low interest 
rate basket, using 6 month blocks. 
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Kendall's T using 6 month blocks for high IR basket 
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Figure 14: Kendall's r for the clayton, frank and gumbel copulas on the high interest rate 
basket, using 6 month blocks. 
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Figure 15: Kendall's r for the clayton, frank and gumbel copulas on the low interest rate 
basket, using 6 month blocks. 
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Figure 16: A"*^'^^^ : 6 month blocks on high interest rate basket. 
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Figure 17: A^'^^^ : 6 month blocks on low interest rate basket. 
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Figure 18: A^^'^^ : 6 month blocks on high interest rate basket. 
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Figure 19: A^^'^^ : 6 month blocks on low interest rate basket. 



39 



j^i,2,3|4 . jgii cjependence for 6 month blocks on high IR basket 

Upper tail dependence 



rVi..^i 



IfXl, ,>l lit Jlni J/^>HlJI III >/vi-^4l j\l lUi\J 



Lower tail dependence 



III In. 



Figure 20: A^^^'^ : 6 month blocks on high interest rate basket. 
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Figure 21: A^^^'^ : 6 month blocks on low interest rate basket. 
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Figure 22: Comparison of Volatility Index (VIX) with upper and lower tail dependence of 
the high interest rate basket. 
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Figure 23: Comparison of Volatility Index (VIX) with upper and lower tail dependence of 
the low interest rate basket. 
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Appendix B 



7.1 Multivariate Clayton Copula 
7.1.1 C^{u) 



7.1.2 ^p'^^ : d-th derivative of the clayton generator 
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7.2 Multivariate Frank Copula 
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7.2.3 Frank Copula Density ( ^ 



d''C 



.du„ 



where 

h^{n) = {l-e-n'-'U{l-e-P-^} 



i=i 



7.3 Multivariate Gumbel Copula 
7.3.1 CG(u) 
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7.4 Multivariate Clayton-Prank-Gumbel Mixture Copula 
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7.4.2 Clayton-Frank-Gumbel Mixture Copula Density 
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